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Constraints on the pairing properties of main sequence 
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ABSTRACT 

Observations of main sequence stars conducted over the last several decades have 
clearly shown that something like 50 per cent of stars of spectral types G and F occur 
in multiple systems. For earlier spectral types, the incidence of multiplicity is even 
higher. Thus, a volume limited sample of white dwarfs should reflect the percentage 
of binarity observed in stars of F to late B spectral types, which are their Main 
Sequence progenitors. However, a study of the local volume limited sample of white 
dwarfs (20 pc from the Sun) conducted by Holberg has shown that a white dwarf has 
a probability of only ~ 32 per cent of occurring in a binary system, in stark contrast 
to the observations of multiplicity of Main Sequence stars. Others studies have also 
led to the same conclusion. 

In this paper, we argue that the "hidden" white dwarfs are either in double white 
dwarf systems or in Sirius-like systems. We also show that the white dwarf progenitors 
of the SDSS white dwarf - M dwarf wide binaries are distributed according to Salpeter's 
IMF. However, they cannot be paired with secondary stars which are also drawn from 
this IMF, since such a pairing would produce a percentage of white dwarf-M dwarfs 
systems that is several times larger than observed. 

Key words: White dwarfs, Binaries: general, stars: formation, stars: low-mass. 
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1 INTRODUCTION 

White dwarf-M dwarf pairs are commonly discovered in any 
surveys of white dwarfs. The observed pairs fall into three 
basic groups. In the first group (group A), the more massive 
star evolves into a white dwarf without ever interacting with 
its companion. These are the so-called wide binaries. In the 
second group (group B), the binary stars have smaller initial 
separations and will eventually evolve through a common en- 
velope (CE) phase as the more massive star becomes a white 
dwarf. These are the post common envelope pre-Cataclysmic 
Variable binaries (Pre CVs). The third group (group C) con- 
sists of pairs that have already evolved through a CE phase 
and are currently seen as close interacting binary systems 
where there is evidence of mass transfer (current or in the 
past). Thus, the M dwarfs in these systems are either filling 
or close to filling their Roche lobes. These are the Cata- 
clysmic Variable (CVs) stars. 

Studies aimed at determining the observed distribution 
M 

/obs(?) of mass ratios q — -—r~ < 1j where M p is the mass 

IVlp 

of the primary star and M s is the mass of its less massive 
companion, are generally limited to establishing the nature 
of the distribution for q ^ 0.1 when the luminosity of the 
secondary star is not swamped by that of the primary. On 
the other hand, insights into the behaviour of /obs(<?) at 



low values of q can be more easily obtained from studies of 
the pairing of white dwarfs with main sequence companions 
belonging to group A and through simple assumptions on 
the initial-final mass relationship for white dwarfs. 

In this paper, we analyse the pairing properties of white 
dwarf - Main Sequence binaries using as constraints the 
white dwarfs - M dwarf sample from th e SDS S belonging 
to group A (wide binaries), the iHolberd ||2009h percentage 
of observed white dwarf binaries in the 20 pc local sample 
and observations of the mass ratio distribution of binary 
Main Sequence stars. We then present some conclusions on 
the pairing properties of Main Sequence low to intermedi- 
ate mass stars, which are the progenitors of the currently 
observed white dwarfs, and link them to possible star for- 
mation scenarios. 



2 THE OBSERVATIONAL BASIS 

2.1 Mass ratio distribution of binary main 
sequence stars 

Early studies of binar y systems point ed to a bi- modal distri- 
bution of / bs(<?) (e.g. lTrimblel (jl97-4T >): a population / bs(g) 
which increases with increasing q and reaches a maximum 
near q = 1 thus showing a preference for stars of similar 
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masses (twins), and a second population /obs(g) which ini- 
tially rises towards lower q, reaches a maximum at q ~ 0.2 
and then decreas es. More recent st u dies u sing F7 to K-type 
primary stars by lHalbwachs et al l (2003) have shown that 
/obs(g) has a very broad peak at q ~ 0.2 — 0.7 (perhaps 
with sub-structure) and a sharp peak for q~0.8 (twins). In- 
terestingly, these authors claim that the f(q) relationship is 
"scale-free" for stars in the three different sub-groups that 
com prised their F7-K star s a mple. 

iKiminki fc Kobulnickvl (|2012) conducted a statistical 
analysis of massive binaries in the Cygnus OB2 associa- 
tion using radial velocity data for 114 B3-05 primary stars. 
They found a mass ratio distribution fobs(q) on q a , with 
a = 0.1 ± 0.5, which is con sistent with their previous work 
( Kobulnic kv fc Frver|l2007|). This is also i n agre ement with 
the observations of iKouwenhoven et al l (|2005l ). who find 
that A and late B-type stars in the Scorpius OB2 association 
have a mass r atio distribution with a = — 0.33, and with 
the studies of IShatskv fc Tokovininl (|2002h . who also sur- 
veyed B-type stars in Scorpius OB2 for binarity and found 
a in the range —0.3 to —0.5. Studies of a sample of mas- 
sive binaries in the Small Ma gellanic Cloud conducted by 
IPinsonneault fc Stanekl (|2006) revealed that their primaries 
appear to have two populations of companions: a "twin" 
population with q > 0.95 comprising 45 per cent of bina- 
ries, and a population with /obs(?) ~ constant (i.e. a = 0) 
comprising 55 per ce nt for their sample. On the other hand, 
ISana fc Evans! (|201ll ) claim that there is no indication for a 
twin population in their Galactic O-star sample and that the 
mass-ratio distr ibution is essenti a lly fla t for 0.2 < q < 1.0. 
The studies by Ragh avan et al l (|2012T ) on companions to 
solar-type stars indicate that the mass ratio distribution is 
flat for ~ 0.2 ^ q ^ 0.95 . 

IWoitas et all (|200lh studied the masses and mass ratios 
of pre-Main Sequence stars. Interestingly, they also found 
that the distribution of mass ratios in T Tauri binaries is 
essentially flat for q > 0.2 and it has no correlation with the 
primary's mass or the stellar separation. They also found 
that there is no significant preference for q ~ 1 (twin com- 
ponent). 

In summary, at the present time it appears that the 
f(q) distribution is either flat or slightly raising toward 
q ~ 0.1 — 0.3 with or without a rise toward q — 1 (twin). In- 
terestingly, these studies appear to agree for primaries rang- 
ing from solar to early O-type stars and even for young pro- 
tostellar associations. However, none of these studies extends 
to low enough mass ratios to investigate the true incidence 
of binaries whose secondaries are M dwarfs. 

2.2 Observations of white dwarf - main sequence 
star binaries 

While an M-dwarf is easily hidden in the glare of an inter- 
mediate or high-mass star primary, such a low luminosity 
star should be more easily seen after its more massive com- 
panion has evolved into a compact star. Thus, insights into 
the behaviour of f(q) at low q can potentially be obtained 
from studies of the pairing of white dwarfs with low-mass 
main sequence stars and by assum ing an initial to final mass 
relationship for white dwarfs fe.g. iFerrario et ail l|2005h ). 

A rather surprising result resides in the ob- 
served incidence of binarity among white dwarfs. 



iFarihi. Becklin fc Zuckermanl (120051 ) found that the 
stellar companion fraction of white dwarfs is 22 per cent, 
uncorrected for bias. Furthermore, they also found that 
most of the stellar compan ions to white dw arfs are low-mass 
M dwarf stars. Similarly, IHolberd (2009) found that only 
32 ± 8 per cent of the local white dwarf population has a 
companion of any type. 

However, it is well known that the progenitors of the 
currently observed white dwarfs, which are Main Sequence 
stars in the mass range ~ 1.2 — 8M©, exhibit a percentage 
of binarity of at least 55 per cent at the low mass limit (F 
type stars), to 60 per cent or more toward the upper end of 
the mass range limit (late B type stars). Thus, there must 
be an additional ~ 30 per cent of as yet undiscovered white 
dwarfs lurking in some kind of binaries. 

Interestingly, if one assumes that the majority of white 
dwarf-M dwarf binaries in the local sample has already been 
detected, since any M dwarf red excess would be clearly visi- 
ble in the spectrum of a near by white dwarf, then according 
to the percentages given bv IHolberd { 2009), the incidence 
of this type of pairing could be as low as ~ 18 per cent, 
consis tent with the findings of lFarihi. Becklin fc Zuckermanl 
(2005). Considering that M-dwarfs are the most numer- 
ous stars in the Galaxy, the fact that they are somewhat 
rarely paired to white dwarfs, suggests that binary forma- 
tion mechanisms tend to exclude them as companions of F 
to late B-type stars. 

Another possibility that cannot be excluded a priori is 
that the local white dwarf binary population is not represen- 
tative of the true Galactic population. If this is the case, our 
calculations, which use as constraint current observations of 
Galactic stellar multiplicity, will simply give predictions that 
would be applicable to a complete white dwarf binary sam- 
ple. 

We show in Fig. fT} the spectral distribution of M 
dwarfs found in wide binaries containin g a white dwarf from 
the sp ectroscopic SDSS Data Release 6 (jRebassa-MansergasI 
l2010h . 

We note that the peak in the M dwarf secondary dis- 
tribution is near M3.5 and that there is a very steep decline 
in the numbe r of M dwarfs of spectral type later than M5. 
Interestingly, IFarihi. Becklin fc Zuckermanl (|2005l ) also find 
that the peak frequency in spectral type occurs around M3.5 
for both field M dwarfs and M dwarf companions to white 
dwarfs (see their Figures 6 and 7). However, relatively to 
the peak, they find that there are ~ 2 — 3 times more L 
dwarfs and ~ 4 — 5 times more M6M9 dwarfs in the field 
than among companions. Thus, despite the excellent sensi- 
tivity in their survey to late-M dwarfs and early-L dwarfs 
very few of these low-mass companions were detected. 

We note that the SDSS white dwarf - M dwarf binary 
data suffer from strong selection effects, since the SDSS is a 
magnitude-limited survey. Thus, one should expect that se- 
lection effects may become dominant at later spectral types 
(M5 to M9) if the white dwarf is more luminous than the 
M dwarf, or at early spectral type (M0 to M2), when the 
fainter white dwarfs may be hidden by the more luminous 
M dwarfs. Therefore, one may expect that even after intro- 
ducing corrections for observational biases, the peak of the 
distribution would remain around M3.5. We shall use this 
peak as a constraint for our studies in section (|4.2[1 . 

An important piece of information that can be extracted 
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Figure 1. Top panel: Spectral clas s of M dwarf companions t o 
white dwarfs in the wide binaries of [Rcbassa-Mans ereas! J201dt) . 
Bottom panel: Mass distribution of white dwarf progenitors after 
a final-to-initial mass relationsh ip has been applied t o the white 
dwarfs in the wide binaries in Rcbassa-Manscrgas (2010) (sec 
text for details). Solid curve: best power law fit with index a = 
-1.95 ±0.14. 



from the data of iRebassa-Mansergas! (|2010l ) concerns the 
distribution of white dwarf masses - and thus of their Main 
Sequence progenitor masses. Since we need an initial-to-final 
mass function to obtain the masses o f the WD progenitors , 
we have used the simple relation of lCatalan et all 12009). 
The histogram of the mass distribution that we have ob- 
tained is shown in the bottom panel of Fig. JT]). This his- 
togram shows that the number of WD progenitors drops 
very quickly with mass. We have found that the best power 
law fit to these data has an index a = —1.95 ± .14, c onsis- 
tent with the observational results of IChabrieij (|2003l ) who 
found a slope a = —2.3 ±0.3 for single stars with M > IMq 
in the galactic disk and young clusters. 

In the next sections, we investigate pairing functions 
of main sequence stars. We will show that given the current 
observational constraints, there should be a sizeable fraction 
of WDs in Sirius-type systems and in double white dwarf 
binaries. 



3 CALCULATIONS 

In order to investigate the properties of white dwarfs in bi- 
naries, a number of systems was generated with primary 
masses distributed according to 

as indicated by the SDSS data (see previous section). The 
justification for this is that the distribution of primary 
masses should exactly reflect the distribution of their masses 



at birth, since in these wide binaries the two stellar com- 
ponents have never come into contact and exchanged mass. 
One should also keep in mind that the mass range of the pri- 
maries is automatically restricted to 1.2~M p ~8Mq since 
the Galactic disk is not old enough for lower mass stars 
to have evolved to the white dwarf stage, while stars with 
M ~ 8Mq undergo supernova explosions and become neu- 
tron stars. 

In our calculations, we have assumed that either (i) the 
masses of the two stars are both drawn from a mass distri- 
bution /(M„) with the constraint M 3 ^ M p (q ^ 1) or (ii) 
the mass of the secondary star depends on the mass of the 
primary, or (iii) the binaries are generated according to a 
Salpeter's like distribution for the primary mass and with 
the secondary mass determined by a generating mass ra- 
tio distribution f(q). Thus we have explored the following 
pairing cases 

(i) A distribution for the mass of the secondaries given by 
dN 



d log M a 



oc MS 



(1) 



where dN is the number of stars per unit volume per 
logarithmic mass inter val d log M s ob served at time t. If 
p = —1.35, we have the ISalpeted l| 19551 ) mass function, 
(ii) A distribution for the mass of the secondaries given by 



dN 



dlogM s 



exp 



Mt 



-(log M s - log Mo) 



2<r 2 



Ms 
M 

1 < 



< 1, 

Ms 

M 



This type o f dist ribution was first 



(2) 

suggested by 



iMiller fc Scale! ll 19791) and was more rec e ntly adopted 
by. e.g.. IChabrieij (|2005l ) and iBochanskil (|2010l ). Here, 
log M c and a are the average mass and standard deviation, 
respec tively, in log M s . In the distribution of Bochanski 
|2010l ) (used in this paper), logM c = -0.745 and a = 0.34. 
(iii) A Gaussian distribution where the mass of the sec- 
ondary is proportional to the mass of the primary star M v : 



dN 



d log M, 



oc exp 



-(logAf B -log(fcAf„)) 



2a 2 



(3) 



where k is the proportionality constant, 
(iv) The mass of the secondary is determined by the mass 
ratio that is drawn from a generating mass ratio distribution 
f(q) given by 

f(q) cc g' 3 < q sC 1 (4) 

We note that if td i s the age of t h e gala ctic disk, assumed 
to be 9.5 Gyr (e.g. lOswalt et all l|l996t )). and T prcw d(M;) 
is the nuclear burning evolutionary time of a star of initial 
mass Mi, then only those stars born at a time ^ r ^ 
td — Tp rcw d can currently be observed as white dwarfs of 
age T w d = td — T prew d — t. In our calculations, we have as- 
sumed that the Galactic star formation rate is is constant 
over the lifetime of the Galactic disc. To obtain an esti- 
mate of r P rewd ( Mi), we have u sed th e average stellar life- 
times given by iRomano et al. I (120051) based on the stellar 
evolution grids of Maeder fc Mevnetl |l989). 

The temperature of the white dwarf can then be cal- 
culated from its age T w d and mass M w d by interpolating 
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Figure 2. From top to bottom: Shaded histograms of d istance, 
effecti ve temperature and mass of the 20 pc sample of iHolberd 
(2009). Dashed histograms: Synthetic white dwarf data. 



the tables for evolutionary sequences of white dwarf atmo- 
sphere s oflHolberg fc Bergeron! ll2006l).lKowalski fc Saumonl 
l|2006t ). lTremblav et all l|201ll ). and iBergeronl ifcOllN H 

Over the evolutionary time of the primary star and sub- 
sequent cooling as a white dwarf, some of the secondaries 
may also have had enough time to evolve and to be ob- 
served as white dwarfs at the present epoch. Thus, in ad- 
dition to systems containing a white dwarf and a Main Se- 
quence star, there will systems comprising of double white 
dwarfs. The mixture of the two will be dictated by the evo- 
lutionary timescales of the secondary star. 

The white dwarf - Main Sequence binaries can then be 
further divided into systems with low luminosity M dwarf 
companions and systems with stars of earlier (K to late B) 
spectral type companions. The latter case yields "Sirius- 
type" systems. If the white dwarf and its main sequence star 
form a close system, then the white dwarf could be hidden in 
the glare of its main sequence companion (just like Sirius B). 
We note that observationally, it is estimated that about 45 
per cent of Sirius-like systems are close or unresolved, with 



1 http: / /www. astro .umontreal . ca/~bergeron/CoolingModels 



the remainder being Close Proper Motion systems (Holberg 
2012, private communication). 

In order to check whether our calculations are consistent 
with the obse rvational char acteristics of the local sample of 
white dwarfs (|Holberdl2009l ), we have plotted in Fig. ([2]) our 
white dwarf synthetic data for distance, effective tempera- 
tu re and mass ov erlapped to the 20 pc white dwarf sample 
of lHolberd l|2009r i. This figure shows that, considering the 
approximations outlined above for the stellar models, our 
synthetic population of non-interacting white dwarfs is con- 
sistent with observational results of isolated white dwarfs. 



4 RESULTS AND DISCUSSION 

4.1 Binary formation mechanisms 

Given the complex nature of star formation, it is difficult to 
predict a priori the expected distribution of mass ratios in 
binary systems, since all stars, including binaries, form in 
clusters with a multitude of physical and dynamical factors 
coming into play. 

Over the years, several routes leading to the formation 
of binaries have been identified. The main ones being (i) 
fission, (ii) capture, (iii) fragmentation. In fission, binaries 
arise from the slow contraction of a rotating gas cloud. Un- 
der the simplified assumption of incompressible, non- viscous 
fluids and of hydrostatic equilibrium, the contraction of a 
rapidly rotating object leads to a dumbbell-type configura- 
tion and to the formation of binaries with mass ratios close 
to unity. This mechanism w a s disp roved by the numerical 
calculations of iDurisen et al.l (1 19861 ) who studied the stabil- 
ity of rotating, compressible gas clouds. They showed that 
the ejection of matter and the formation of a bar with trail- 
ing spiral structures can efficiently redistribute the angular 
momentum on a dynamical time scale avoiding the forma- 
tion of a binary system. 

The breakup of a fast rotating gas cloud while it is 
collapsing is referred to as "pro mpt fragment ation" . This 
mechanism was first proposed by iHovla l|l95&t ). Because of 
low compressional heating, the gas can cool radiatively and 
the contracting gas cloud is initially isothermal. As the cloud 
collapses, the density increases and the Jeans mass decreases 
causing fragmentation to take place. However, as contrac- 
tion continues, the gas cloud's opacity increases causing its 
core to heat up. Fragmentation ends once a fragment's tem- 
perature, and thus its Jeans' mass, starts to increase. The 
critical value of opacity that halts contraction yields a min- 
imum Jeans' mass value of about O.O1M0. Early computa- 
tional work on the forma tion of binary systems via fr agmen- 
tation w as cond u cted b y lBoss fc Bodenheimerl ljl979l ). More 
recently, I Clarke! (| 19961 ) studied the scale-free fragmentation 
scenario and its observational implications. 

Fragmentation of a protostellar disc, also called "disc 
fragmentation", occurs once the collapse of the cloud 
is over and p r otoste llar objects have already formed. 
iBonnell fc Bate! (1 19941 ) proposed that fragmentation can oc- 
cur via two mechanisms. One involves rotational instabili- 
ties in a protostellar core. The other invokes the formation 
of a rapidly rotating core which is unstable to axisymmetric 
perturbations. This core bounces into a ring which quickly 
fragments into several components. 
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The capture mechanism implies that single stars have 
already formed and are still in a (dense) stellar cluster. Bina- 
ries are then created when one star captur e s ano ther ("dy- 
namical capture", see iMcDonald fc Clarke] (|l993h >. The re- 
quired excess kinetic energy loss could occur via tidal dis- 
sipation, if the stars get sufficiently close to each other, or 
through dynamical encounter and transfer of energy to a 
third star. If capture occurs early on in the formation process 
and either one or both young stars are still surrounded by 
protostellar discs, then energy dissipation would effectively 
occur through discs-stars interaction ("star-disc capture", 
IMcDonald fc Clarke! (|l995l )V 

Clearly, most of the above mechanisms are expected to 
play a role in a star forming region, although it is gener- 
ally expected that cloud fragmentation may play a dom- 
inant role. This is certainly the picture that is emerging 
from the hydrodynamic simulatio ns of s tar cl uster forma- 
tion, such as those conducted by iBatei (|2012T ). which can 
resolve masses dow n to a few Jupi ter masses. Their calcu- 
lations produce the lChabrierl (|2005r i observed IMF and find 
that stellar multiplicity increases with the mass of the pri- 
mary, in agreement with observations. Furthermore, they 
find a mass ratio distribution for solar-type and M-dwarf 
binaries which is roughly flat, again consistent with obser- 
vations. 

It is interesting to note that the flat fob B {q) distribution 
observed in protostellar objects seem to support the view 
that in most multiple systems in T Tauri associations the 
components' masses are mainly determined by the cloud's 
fragmentation process and not by the subsequent disc accre- 
tion processes (|Woitas et alj|200l] ). This finding is corrobo- 
rated by the numerical simulations of lBatei l|l997h who have 
shown that in the disc accretion stage following cloud frag- 
mentation, the system's mass ratio would tend to approach 
u nity, which is i n disa greement with the observational work 
of lWoitas~e t al. (200l]). This result is very significant, since it 
indicates that binary characteristics are already determined 
as early as 1 million years after cloud fragmentation has 
begun. 

We stress that what is observed (and usually fitted 
with power laws), is the "specific mass ratio" distribution 
fobs(q) which is obtained for a sample of stars with primaries 
in a g iven mass range {M v \ and M P 2, iKouwenhoven et al.l 
(2009)). This is certainly the case if one studies the pairing 
properties of binaries whose primary stars are white dwarfs, 
since the progenitors of white dwarfs are main sequence stars 
with masses 1.2$M P S 8M . The observed mass ratio dis- 
tribution is also usually confined to g~0.1 because observa- 
tions are often incomplete below this value. 

In what follows, we will consider the pairing possibilities 
listed in section ((3]) to see whether the existing observational 
constraints, comprising the SDSS white dwarf - M dwarf 
wide binaries, the local white dwarf binary sample, and the 
flatness of the observed mass ratio distribution, favour any 
of the aforementioned routes. 

4.2 Results of calculations and comparison with 
observations 

In the stellar capture scenario, single stars in dense, young 
clusters form binaries via "random pairing". This can be 
simulated if we draw the masses of both stars from an IMF 



distribution of the type given by ((2]). The result of such pair- 
ing yields ~ 85 per cent of white dwarf - M dwa rf systems, in 
disagreement with the study of iHolberg] (|2009T l which gives 
an upper limit of ~ 40 per cent (see section (|2.2[) '). Such a 
distribution would also give a peak in the mass ratio distri- 
bution near q = 0.2 and a very steep decline toward 9=1, 
which is in disagreement with the observed, flat behaviour 
for 0.1~g~l discussed in section [2~T1 We show the results 
in Fig. In this figure we also show the predicted distribu- 
tion of secondary star masses (centre panel) and M-dwarfs 
(right-end panel) . We note that the peak of the M dwarf dis- 
tribution is near spectral class M6, while the SDSS peak of 
the M dwarf distribution is near M3-M4. We also note that if 
both primaries and their companions are selected randomly 
from the IMF, then the outcome would be a small percent- 
age of Sirius-type systems (~ 15 per cent) and double white 
dwarfs (~ 2 per cent). 

We have also explored a distribution as given in (0. 
Here, the type of companion depends on the mass of the 
primary. The behaviour of /obs(g) and M s can exhibit peaks 
at different locations and widths depending on the values 
assigned to k and a. If we set a = 0.4 and k = 0.5 we obtain 
a low percentage (~ 20 per cent) of white dwarfs - M dwarf 
systems and a high percentage of Sirius-type and double 
white dwarf systems of 57 and 23 per cent respectively. Fig. 
13] shows the results for this set of values. We also note that 
the peak of the M dwarf distribution for the companions is 
predicted to occur around M3-M4, as observed. However the 
slope of the mass ratio distribution for g~ 0.1 — 0.2 is again 
too steep to fit the observational evidence pointing to a flat 

/ob S (g)- 

In the pairing cases considered above, the masses of the 
two stars have been drawn from some mass distributions 
which yields binaries with a certain observed mass ratio dis- 
tribution. What we are going to look at now, is the situation 
where the physical processes leading to binary star forma- 
tion establish the distribution of the primary mass and of 
the mass ratio distribution /(g) (the "generating mass ratio" 
distribution) . This approach on generating binaries could be 
interpreted in the framework of t he scale-free f ragmentation 
scenario for binary formation of I Clarke] (Il996l ) who consid- 
ered a collapsing and breaking up gas cloud whose final frag- 
ments are clumps of mass of M c i ump which are distributed 
according to a mass function -F c i ump (M c i ump ). In her model, 
each clump then divides into two pieces whose mass ratio 
distribution is /(g). That is, the fraction of clumps with 
mass ratio in the range q to q + dq is f(q)dq. In her stud- 
ies, the function /(g) is independent of M c i ump , which is the 
necessary assumption for scale-free cloud fragmentation. 

In our calculations, we have considered a uniform mass 
ratio distribution /(g) and a power-law for the primary mass 
distribution f(M p ), as indicated by the SDSS data. It is im- 
portant to note that the stellar mass limits play a crucial role 
in the determination of the observed mass ratio distribution 
/obs(g). In the studies of white dwarf binaries, however, the 
spectral type of the white dwarf's progenitor is restricted to 
be between F and late B which gives rise to the flat /obs(g) 
shown in Fig. ((5j). Such a pairing yields a small percentage of 
white dwarfs with M dwarf companions (~ 18 per cent) and 
a large percentage of Sirius-type systems (~ 47) and double 
white dwarfs (~ 35). The peak of the M dwarf distribution 
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Figure 3. From left to right: Mass ratio distribution, companion mass distribution f c {M a ) and spectral class distribution of M dwarf 
companions for a pairing where the masses of the components are drawn from the IMF of the type The panel on the right shows 
the normalised SDSS M dwarf distribution of wide white dwarf-M dwarf binaries (dashed line), while the solid line shows the model 
prediction for the distribution of M dwarfs. 




Figure 4. From left to right: Mass ratio distribution, companion mass distribution fc(M s ) and spectral class distribution of M dwarf 
companions for a pairing of the type $3^ with u = 0.4 and k = 0.5. The panel on the right shows the normalised SDSS M dwarf 
distribution of wide white dwarf-M dwarf binaries (dashed line), while the solid line shows the model prediction for the distribution of 
M dwarfs. 



also falls near M4, which is in general agreement with the 
SDSS observations. 

Finally, it may be of interest to address the fact that 
some observations show the existence of a "twin peak", 
caused by stars of similar masses pairing up. There have 
been some doubts on whether this peak is real, since bi- 
naries whose components have similar masses tend to be 
brighter than those with fainter companions, thus result- 
ing into a poss i ble ov ersamplirig of these sys tems (e.g. 
lHalbwachs et alj <l2003h'). On the other hand, iTokovininl 
(2000.;) and lHalbwachs et all (|2003l ) also report that the fre- 
quency of twins does seem to be higher at short orbital pe- 
riods (< 40 — 50 days). They also note that the mass ratios 
of their visual binary sample with wide separat ions does not 
exhibit a peak around q = 1 (|Halbwachslll98 j ) . 

The "twin peak" seems also to be c o mmon among mas- 
sive primaries. iPinsonneault fc Stanekl |2006l ) report that 
systems with M2 > 0.95Mi comprise 45 per cent of their 



population of massive binaries while the rest exhibits a flat 
mass ratio distribution. 

Thus, it appears that the stellar components of binaries 
with short periods and some of the more massive binaries 
"prefer" similar mass companions. This implies the coexis- 
tence of different binary formation mechanisms. 

One possibility is that in some cases, dictated by some 
as yet unknown initial conditions, the scale-free fragmenta- 
tion of the cloud into clumps is then followed by vigorous 
disc fragmentation (see section l4"T1) . As theoret ically demon - 
strated by the hydro-dynamical calculations of lBatd |l997l ). 
such disc fragmentation would create systems whose mass 
ratios approach unity. This behaviour can be modelled by 
setting k = 1.0 in equation ([3|. In this picture, scale-free 
fragmentation would then be followed by a formation pro- 
cess that is dependent on the mass of the primary star. 

In Table 14.21 we summarise our results in terms of the 
various systems that are generated via the pairing mecha- 
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Figure 5. From left to right: Mass ratio distribution, companion mass distribution f c {M a ) and spectral class distribution of M dwarf 
companions for a pairing of the type {JJl with = 0. The panel on the right shows the normalised SDSS M dwarf distribution of wide 
white dwarf-M dwarf binaries (dashed line), while the solid line shows the model prediction for the distribution of M dwarfs. 



nisms considered in this paper. This table gives a prediction 
of the relative percentages of systems that one should find in 
a complete volume limited sample of white dwarf binaries. 



5 CONCLUSIONS 

In this paper, we have seen that it is possible to gain some 
insights in the behaviour of the mass ratio distribution at 
low values of q by studying binary systems comprising of 
white dwarf primaries and main sequence companions. This 
has allowed us, via an initial-final mass relationship for white 
dwarfs, to draw some conclusions on the pairing properties 
of Main Sequence of spectral type F to late B, which are 
the progenitors of the currently observed white dwarfs. The 
constraints that we have used are (i) the white d warfs - M 
dwarf wide binary sample from the SDSS (ii) the iHolberd 
(2009) percentage of observed white dwarf binaries in the 
20 pc local sample and (iii) the fact that the mass ratio 
distribution for Galactic binaries, ranging from G to early 
B-type stars and even for young protostellar associations, is 
either flat or slightly raising toward q ~ 0.1 — 0.3. The aim 
of this paper was to investigate whether it is possible, on the 
basis of the data currently at hand, to explore and possibly 
favour some of the proposed routes for star formation. 

A prediction of our studies is the existence of a large 
fraction of "hidden" white dwarfs in binaries consisting of 
either double white dwarfs or close (unresolved) Sirius-like 
systems. 

In this context, we would also like to note that using the 
ROS AT Wide Field Camera survey of th e extreme ultravi- 
olet, |Burleigh^^arst^^^£olber3 (|l998l ) have revealed the 
existence of a previously unidentified sample of hot white 
dwarfs in unresolved, detached binary systems. These stars 
are invisible at optical wavelengths due to the proximity 
of their much more luminous companions (spectral type K 
or earlier). However, for companions of spectral type A5 or 
later, the white dwarfs are easily visible at far-ultraviolet 
wavelengths, and can be identified in UV spectra. In total, 
16 such systems have been discovered in this way through 
ROSAT, EUVE and IUE observations. Further observations 



of this kind should reveal how common these systems really 
are. 

Clearly, also a much larger sample of white dwarf - M 
dwarf companions is needed to further explore the incidence 
of such systems. At the moment, the number is still too 
small and uncorrected for observational biases. A clean, en- 
larged sample of this type is crucial to shed more light on 
how binaries are formed and thus on possible star formation 
mechanisms. 
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